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Abstract 

A non-subtractive recipe of Casimir energy renormalization efficient in the pres- 
ence of logarithmically divergent terms is proposed. It is demonstrated that it can 
be applied even then, when energy levels can be obtained only numerically and 
neither their asymptotical behavior, nor the analytic form of spectral equation is 
known. The results of calculations performed with this method are compared to 
those obtained by means of explicit subtraction of divergent terms from energy. 

1 Introduction 

Ever since Casimir [1] has obtained corrections to the energy of a macroscopic 
system due to vacuum fluctuations of quantized electromagnetic field in 1948 
this effect has been intensively studied both from theoretical and experimental 
points of view. Nevertheless the calculation of Casimir energy except for the 
most simple problems involving free fields inside cavities with flat boundaries 
is quite non-trivial yet. To realize this, recall a great number of papers devoted 
to a free field in the interior of a sphere [2-11]. Despite the fact that in this 
case one can explicitly write out spectral equations, the first analytical results 
have been obtained only in [12] for massive scalar field and in [13] for fermions. 

It should be stressed that the knowledge of analytical form of spectral 
equation has been crucial in [12] since it makes possible the transition from the 
sums containing the unknown energy levels to the integrals with the explicit 
integrands [14]. The goal of this paper is to demonstrate a method which 
can be applied to numerical calculation of Casimir energy in cases when these 
requirements are not met. Moreover, we are not going to use a rather standard 
trick [7], [15], which lets one overcome problems arising due to the presence of 
logarithmic divergency in Casimir energy of free massless fields inside spherical 
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shells. Note that logarithmic divergency appears as a consequence of a curved 
surface bounding the shell and makes the energy renormalization ambiguous. 
The main idea of the trick is to consider the "inner" and "extra" problems 
together since their logarithmic divergencies cancel each other. 

The ambiguity of Casimir energy renormalization in the presence of log- 
arithmic divergency is quite obvious. Indeed, in case of massless fields the 
energy of the system can be characterized by a single dimensional parameter 
L which is the linear size of the system. The regularization parameter a can 
be also chosen to have a dimension of length. In the absence of logarithmic 
divergency the "minimal subtraction" of singular terms is not only natural 
but also well grounded. Indeed, any term proportional to a~ s (s > 0) is obvi- 
ously proportional to L s_1 , i.e. to the non-negative power of L. This makes it 
possible to normalize the final result at L = oo, where Casimir energy should 
become zero, and subtract all singular terms at the same time. After such 
subtraction the only remaining term in the limit a — > 0, which reads ca°/L, 
provides the final result. 

In the presence of logarithmic divergency the subtraction becomes am- 
biguous, since in order to renormalize the term ca° log(a / L) / L , one should 
subtract ca°\og(da/L)/L, where d is an arbitrary constant, which cannot be 
determined from the normalizing condition at L — > oo. 



2 Massive scalar field in ID 

To illustrate the main idea of the proposed recipe let's consider Casimir energy 
with the logarithmic divergency in the most simple case, i.e. Casimir energy of 
massive scalar field on interval of length L with Dirichlet boudary conditions 
at the ends of the interval. It reads: 

£cas = - ^ = - \](™/L) 2 + m 2 (1) 

A n=l Z n=l 

Note that we aren't going to subtract the Minkowski vacuum contribution 
from (1) as proposed in "standard" approaches. Instead we will directly pass 
to renormalizaton. 

The regularization of (1) requires the introduction of the parameter a which 
has a dimension of length and stands in the argument of the cut-off function 
F(au n ): 

1 oo 

= ~ E UnF{au n ) (2) 

Z n=l 

Trivial considerations based on dimensional analysis lead to the following ex- 
pression for the regularized Casimir energy 

L L° 1 
Zcas - c -2^ + c_i— r + c Q - + c x m 2 L log(a/L) H (3) 

It can be easily verified that for various cut-off functions such as F(x) = 
exp(— x), F(x) = exp(— x 2 ), F(x) = exp(— x 3 ), . . ., F(x) = exp(— x 6 ), F(x) = 
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exp(— 2 cosh(x) +2), . . ., identical c\ are obtained, while Cq are different. The 
identity of c\ for different F(x) can be demonstrated with the following esti- 
mation for the sum giving rise to the logarithmic divergency: 

2£^ fK) ~ ^logJV- — logfAsL/a) = 

— log(L/a) + — log(Ax) , (4) 

where N ~ AxL/a and Ax is a cut-off interval of F(x). 

Since any subtraction in the presence of logarithmic divergency is am- 
biguous this procedure should be excluded from consideration along with the 
logarithmic divergency itself. To realize that let's calculate d\£^ s : 

00^ + 0^ + ... (5) 

The obtained expression is regular in the limit a — > 0, so no subtraction 
is required. The knowledge of the function df£<£2, lets one reconstruct the 
required £^ s unambiguously, since the initial conditions at L — > oo are well- 
known: both £^\ and 0l£^1 should become zero. Note that while (5) doesn't 
describe the asymptotical behavior of d^£^l at L — > oo, it demonstrates 
the disappearence of all singular terms in it. Moreover the following integral 
approximation shows that d^j^l vanishes for L — > oo: 



1 00 1 00 m 

£cl = 9 E ^ n F(auj n ) = jJ2 u n F(au n ) - —F(am) « 

Z ji=1 4 n=-oo 4 

~ - / dxih I tx)\J x 1 + m 2 F(aVx 2 + m 2 ) — -F(am) (6) 
2 J— oo 4 



3 Method of calculation in general case 

Let's generalize the proposed method in such a way that it doesn't require 
the analytical expression for energy levels. Suppose one has a set of energy 
levels of some spectrum u n and the corresponding Casimir energy contains the 
logarithmic divergency. First of all it turns out to be possible to modify the 
initial expression for the Casimir energy by introduction of some parameter /i 
in such a way that 

1 00 / / 

~ E ^l-^F{a^l-^) (7) 

Z n=l 

doesn't contain the logarithmic divergency. For the massive scalar field on 
an interval \x is obviously equal to the mass of the field. In less trivial three- 
dimensional cases with spherical symmetry fi is some parameter having a 
dimension of mass which characterizes the total coefficient by the logarithmic 
divergency with all values of angular momentum taken into account. 

The next step is to introduce an "additional mass" of the field M. and 
study the Casimir as a function of it in the range from M. = to M. = oo. In 
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fact it's helpful to introduce another parameter M: M 2 = M 2 + /i 2 and study 
the modified Casimir energy 



1 00 / 

£UM) = - £ \R + M 2 -^ 2 (8) 

A n=l 

as a function of M in the range from M = 00 to M = /i. To realize this one 
should calculate numerically the following quantity in the specified range of 
M: 



d 2 M (££ S (M)/M 



d 2 M 




(9) 



Note that in contrast to our previous considerations we have substituted di- 
mensionless quantity 



f^»+ M2 -" 2 =Vt + 1 

to the argument of so that ct should be also taken dimensionless. 

An alternative interpretation of (9) follows from the observation that M 
acts as an effective length L in the expression for Casimir energy. Indeed, 
(9) can be obtained from the initial expression as a result of the folowing 
transformation of the spectrum. At the first step u n is transformed to u' n = 
\Joj~ 2 , — fi 2 / fi, which is dimensionless spectrum with the subtracted effective 
mass. After that the scale transformation of the system x — > x(M/fj) and 
uo' n — > id" = u' n /(M/fi) is performed. In the end the unit mass is "added" to 
the obtained spectrum: 



The limit M — > 00 obviously corresponds to the infinite size of the system, 
while for M — /i one obtains the initial spectrum divided by //. 

It's easy to see that all divergent terms in (9) vanish. In the limit n — > 00 
one can make use of the following expansion 



(11) 



The first term in this expansion gives rise to the sum which is free of logarith- 
mic divergency due to the definition of /i. Other divergencies are proportional 
to (M/a) 2 /M and (M/aY/M and vanish when the second-order derivative 
is calculated. The second term leads to the logarithmic divergency with the 
coefficient proportional to M 1 by it, which also vanishes. 

As a result one has (9) regular for a — > and the natural normalizing 
condition £ cas (M — > 00) = 0. The latter can be understood from two different 
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points of view. On one hand the quantized field with the infinitely large mass 
should have zero Casimir energy. On the other hand the Casimir energy in 
the limit of the infinite size of the system should become zero. Whichever 
interpretation is chosen, the obtained results let one reconstruct the required 
S cas (M = //) which corresponds to the initial spectrum. 

Note, that principally one could consider £&£(M) instead of £M(M)/M. 
However that would increase the order of derivative required to exclude all 
divergent terms by one what is undesirable from the practical point of view. 

The proposed method turns out to be efficient not only in the most trivial 
one-dimensional cases but also in more realistic three-dimensional ones. How- 
ever to employ it in three-dimensional case one should inevitably calculate the 
fourth-order derivative of the Casimir energy (9) since the main singular term, 
which is proportional to the volume of the system, reads C-^L 3 /a 4 . It should 
be also noted that in this case the calculations of the sums become more so- 
phisticated since the final value of a typical sum is about 40 orders lower than 
intermediate values obtained during its calculation and extra floating-point 
precision is required. 

4 Numerical results 

For scalar field on an interval [0; L] with L — 1 spectrum reads 



The result of the straightforward application of our method with vari- 
ous cut-off functions such as F(x) = exp(— x), F(x) = exp(— x 2 ), F(x) = 
exp(— x 3 ), . . ., F(x) = exp(— x 6 ), F(x) = exp(— 2 cosh(x) + 2) is presented on 
Fig.l. It has been shown that for each of these functions the same result is 
obtained and what's more the precision of coincidence depends only on the 
number of energy levels taken into account and the number of right digits used 
in the realization of floating-point arithmetics as well. 

As to dependence of the Casimir energy on the mass of the field some 
important aspects should be stressed. First of all in the limit m — > a 
well-known result for the massless scalar field is obtained. In the range of 
large values of m Casimir energy decreases exponentially as e~ 2mL what could 
be expected from qualitative considerations. The results obtained with our 
method in this trivial case are in agreement with those obtained using the 
traditional subtractive technique. 

To demonstrate how the method can be employed in less trivial cases we 
consider the massless scalar field inside of a spherical shell of radius R — 1. 
In this case the same set of cut-off functions has been used. The values of 
an effective mass fi = 0.1377 obtained with each of these functions coincide 
up to the first four digits. Consequently the precision of the obtained E cas = 
3.790 • 10~ 3 has the same order, what corresponds to about 200 s-levels taken 
into account during the calculations. The number of energy levels taken into 
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account is directly affected by the range which the regularization parameter a 
used in the calculations belongs to. Therefore one can control precision of the 
final result simply changing the range of employed values of the regularization 
parameter. 

Note that in the framework of this approach we have obtained not only 
Casimir energy of the massless scalar field (corresponding to M. = 0) but also 
Casimir energy for all possible values of mass in the range from zero to the 
"effective" infinity. The dependence of the Casimir energy of the scalar field 
inside the sphere on the mass of the field is presented on Fig. 2. 

It should be stressed that on the contrary to the results obtained with the 
traditional subtractive technique in [12] our result doesn't contain logarith- 
mical singularity at M. = what seems more reasonable from physical point 
of view. The most likely explanation of this difference is that the subtractive 
procedure contains some arbitrariness. As a result some function which has 
regular behavior at M. — > oo but is singular at M. — > could be subtracted 
from the final result. 

As has been pointed out in [12,15] there is no argument at present which 
can remove this arbitrariness in case of a massless scalar field inside of a 
sphere. Therefore Casimir effect in the whole space with Dirichlet boundary 
conditions on the sphere is usually considered. It seems reasonable to calculate 
the Casimir energy in the last case employing our method. 

In fact there are two ways to proceed to take exterior into account. The first 
one deals with the continuous spectrum and requires that all the regularized 
sums be replaced with the appropriate integrals containing the energy levels 
density in the integrand. The second way lets one work with discrete spectrum 
all the time. To realize that one should place the initial spherical shell into 
another sphere with the radius R ou t = kRi n where k > 1 and calculate the 
Casimir energy for the system bounded by the outer sphere taking into account 
boundary conditions on both of the spheres. For each finite k the spectrum 
is discrete and the developed technique can be applied without modification. 
The required result can be achieved in the limit k — > oo. In practice it turns 
out that for k > ko, where ko is finite and depends on the required precision, 
the result doesn't depend on k. It turned out that in the considered case in 
order to get 4 right digits in the final result k ~ 10 is quite enough. 

The final result of the calculations is presented on Fig.3. Note that while 
the qualitative behavior of Casimir energy is the same as that obtained with 
methods employing explicit subtraction [15], there is no absolute coincidence. 
For example, for the massless field the result obtained with our recipe is 
£ cas (M = 0) = 0.0039 while direct subtraction leads to £ cas (M = 0) = 0.0028. 

5 Conclusion 

To summarize, an efficient technique for numerical calculation of Casimir en- 
ergy in the presence of logarithmical divergencies has been developed. The 
advantages of the proposed method are its ideological triviality and universal- 
ity which let one apply it to a wide range of problems in which numerical values 
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for energy levels can be obtained. The results of its application to a number of 
problems appear to be reasonable, especially in case of a curved boundary. As 
to disadvantages they are purely technical: one should employ floating-point 
arithmetics with extra precision to carry out calculations in realistic cases. 
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Fig.l. The Casimir energy of the massive scalar field on a unit interval as a function of 
the mass of the field. 
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Fig. 2. The Casimir energy of the scalar field inside of a spherical shell of radius 1 obeying 
Dirichlet boundary conditions as a function of the mass of the field. 
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Fig.3. The Casimir energy of a massive scalar field in the whole space with Dirichlet boundary 
condition on a sphere of radius 1 as a function of the mass of the field. 
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